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Abstract. We study a family of the Laurent biorthogonal polynomials arising from the Her- 
mite continued fraction for a ratio of two complete elliptic integrals. Recurrence coefficients, 
explicit expression and the weight function for these polynomials are obtained. We construct 
also a new explicit example of the Szcgo polynomials orthogonal on the unit circle. Relations 
| with associated Legendrc polynomials are considered. 
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1 Introduction 

We start with well known identity (see e.g. |24|) 



d f 1 . - = n(n - 1) sn n " 2 u - n 2 (l + k 2 ) sn n u + n(n + l)k 2 sn n+2 u (1.1) 
du z 

for the Jacobi elliptic function sn?j depending on an argument u and a modulus k. This (and 
similar 11 formulas obtained if one replaces sn with other elliptic functions) identity goes back 
to Jacobi and usually is exploited in order to establish recurrence relations for elliptic integrals. 
Indeed, introduce the following elliptic integrals 

rx L,2n-i2n Jj. 

2n „ 2n ' tv L UL 







J n (x;k) = / k zn sn' n udu= / , (1.2) 

Jo V(l-* 2 )(1-* 2 * 2 ) 



where snn = t and x = sni>. Then from (jl.lj) we obtain 

(2n - 1) J n (x; k) - (2n - 2)(A; 2 + 1) J„_i(x; k) + A: 2 (2n - 3) J n _ 2 (x; k) 



*This paper is a contribution to the Vadim Kuznetsov Memorial Issue "Integrable Systems and Related Topics". 
The full collection is available at |http://www.emis.de/journals/SIGMA/kuznetsov.html| 
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= k 2n - 2 x 2n ~ 3 ^(l - x 2 )(l - k 2 x 2 ). (1.3) 
From this formula we can express J n +i(x; k) in the form [12] 



J n+l {x-k) = Q(x; n)y/(l - x 2 ){\ - k 2 x 2 ) + A n Ji{x] k) - B n J {x; k), (1.4) 

where Q(x;n) is a polynomial in x (depending on n) and the coefficients A n , B n satisfy the 
same recurrence relations as J n+ i(l; k), i.e. 

(2n + l)A n - 2n(l + k 2 )A n ^ + (2n - l)fe 2 A n _ 2 = 0, 

(2n + l)B n - 2n(l + k 2 )B n ^ + (2n - l)fc 2 fl„_ 2 = (1.5) 

with obvious initial conditions 

A_i = 0, B-i = -l, A = l, S = 0. (1.6) 

Note that from recurrence relations fll.5|) and initial conditions (jl.6p it follows that ^4 n (A;) is 
a polynomial in A; 2 of degree n and B n (k) is a polynomial in A; 2 of degree n but having common 
factor k 2 , i.e. B n = k 2 V n -i(k 2 ), where V n -\(z) is a polynomial of degree n — 1 in z for any 
n=l,2,.... 

There is an elementary Wronskian-type identity following directly from fll.5|) |12j : 

Zj,2n 

B n A n -i — A n B n _i = - — — . (1.7) 
in + 1 

Formula (|1.4|) allows to reduce calculation of any (incomplete) elliptic integrals of the form 
Pit) dt 



X 



y/(l-t 2 )(l-k 2 t 2 ) 

(where P(t) is a polynomial) to standard elliptic integrals of the first and second kind Jq(x; k) 
and J\(x;k). This result is well known since Jacobi. 

In what follows we will denote Jo(l;fc) = K(k) and J\(l;k) = J(k). Note that K(k) is the 
standard complete elliptic integral of the first kind [24] and 

dE 

J{k) = K(k) - E(k) = -k=r, 

dk 



where 



1 'l-ft^ 1 ' 1 



is the complete elliptic integral of the second kind. 

We note also an important relation with hypergeometric functions 

K{k) = | 2 F!(l/2, 1/2; 1; k 2 ), E{k) = | 2 F 1 (-l/2, 1/2; 1; k 2 ), 

J(k) = ^f 2 F 1 (3/2,l/2;2;k 2 ). (1.8) 

Hermite in his famous "Cours d'analyse" [12] (see also \10\ lllj) derived a continued fraction 
connected with a ratio of two complete elliptic integrals: 

J(k) k 2 



K(k) n/< _ 9k 2 



;i.9) 



2(1 + r 



4(1 + t) -6(TTA?)" 
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where 



y/(l-X 2 )(l-k 2 X 2 ) JO y/(l-X 2 )(l-k 2 X 2 ) 

The Hermite continued fraction (|1.9|) follows directly from (|1.3|) . Note that the continued frac- 
tion (|1.9p belongs to the class of the so-called T-continued fractions (with respect to the variab- 
le k 2 ) introduced and studied by Thron (see [Uj for details). Thus perhaps Hermite was the 
first to introduce an explicit example of the T-continued fraction. In what follows we will see 
that this example gives rise to a class of polynomials which are biorthogonal on the unit circle. 
Rewrite relation (|1.4p for x = 1 in the form 

J(k)/K(k)-B n /A n = ^^- (1.10) 

from which Hermite concluded that the rational function B n /A n is an approximate expression 
for the ratio J(k)/K(k) to within terms of degree n + 1 in k 2 . 

Hermite also noted that the coefficients A n , B n appeared as power coefficients in the following 
Taylor expansions: 

J ° (X) -A x + A l x 3 + --- + A n x 2n+1 + ---, 



v/(l-x 2 )(l-fc 2 x 2 ) 

= B x + Bix* H + B n x zn+1 H . (1.11) 



^(1 ~X 2 )(1 - k 2 X 2 ] 

Thus formulas (jl.lip can be considered as generating functions for A n , B n . 



2 Laurent biorthogonal polynomials 

Introduce the new variable z = k 2 and define the polynomials P n {z) = A n /^ n of degree n in z, 
where 

^ 2n + 3 

?o = I- 



^n+i 2n + 2' 
We have 

n\ 

where (a) n = a(a + 1) • ■ ■ (a + n — 1) is the standard Pochhammer symbol (shifted factorial). 

Then it is seen that P n {z) = z n + 0(z n_1 ), i.e. P n {z) are monic polynomials. Moreover 
from (jl.5p it follows that P n {z) satisfy the 3-term recurrence relation 

P n+1 (z) + d n P n (z) = z(P n (z) + W„_i(s)), (2.2) 

where 

fn + 1/2) 2 

d n = -l, n = 0,l,..., b n = - \ /-, , n = l,2,... (2.3) 

n[n + 1) 

with initial conditions 



P_i(z) = 0, fb(z) = l. 



(2.4) 
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Define also the polynomials 
P^ 1 (z) = 2Bj(zt n ). 

It is then easily verified that the polynomials Pn\z) are again n-th degree monic polynomials 
in z satisfying the recurrence relation 

Pi%(z) + d n+l P^(z) = z(P£\z) + 6 n+1 P« 1 (^)). (2.5) 

The polynomials P n (z) are Laurent biorthogonal polynomials (LBP) [7] and the polynomi- 
als Pn\z) are the corresponding associated LBP. The recurrence coefficients b n , d n completely 
characterize LBP. The nondegeneracy condition b n d n ^ 0, n = 1, 2, . . . [7J [27] obviously holds 
in our case. 

It is well known that LBP possess the biorthogonality property [7J. This means that there 
exists a family of other LBP P n (z) and a linear functional a such that 

(a, P n {z)P m {l/z)) = h m 5 nm , (2.6) 

where the normalization constants h n are expressed as (7j 

k=i k 

The linear functional a is defined on the space of all monomials z s with both positive and 
negative values of s: 

c s = {a,z s ), s = 0,±l,±2,..., (2.7) 

where c s is a sequence of moments (this sequence is infinite in both directions). 

The biorthogonality condition ()2.6[) is equivalent (under the nondegeneracy condition h n ^ 0) 
to the orthogonality relations 

{a,P n (z)z- j )=0, j = 0,l,...,n-l. (2.8) 

Note that in our case from (j 2 . 2 j) and (|2.3p it follows that 

P n (0) = l, n = 0,1,2,.... (2.9) 

Introduce the reciprocal LBP by the formula 



It appears that P^iz) are again LBP with the recurrence coefficients [23] 

< = ^- (2.H) 



The moments c* of the reciprocal polynomials are expressed as [23] 

< = (2-12) 
ci 

In our case fj2.3j) we have 6* = 6 n , d* n = d n . Hence the reciprocal polynomials coincide with 
the initial ones: 

z n P n (l/z) = P n {z). (2.13) 
Moreover in our case c\ = do = — 1 and thus c* = —c\- n . 
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The polynomials P n {z) are the biorthogonal partners with respect to the polynomials P n {z). 
Their moment sequence c n is obtained from the initial moment sequence by reflection: 

c n = c_ n , n = 0,±l,±2,.... (2.14) 

Note that under the assumption cq = 1 (standard normalization condition) we have Co = 1 as 
well. 

For the biorthogonal partners there is an explicit expression [23] 

Pn[z) - «j • (2 - 15) 

In our case (taking into account properties (|2.9p and (j2. 13[) ) it is seen that 

p n{z) = Pn + l(z)-P n (z) _ (2 lg) 

Formula f|2. 16j) admits another interpretation if one introduces the Christoffel transform (CT) 
of LBP. Recall gTJ [23J [22] that the CT for LBP is defined as 

P iC) {z) = P n + ^)-UnPn(z) ^ ^ = Q> ^ U n = ^^, (2.17) 

(C) 

where /i is an arbitrary parameter such that P n {n) ^ 0, n = 1, 2, . . . . The polynomials Pn (z) 
are again monic LBP with the transformed recurrence coefficients 

On. + «n + l^n 



(2.19) 



The moments ch corresponding to CT are expressed as [27] 

(C) _ °n+l ~ fJ'Cn 
Cl - /J, 

There is a special case of the CT when fi = 0. In this case, we have 

p (C) {z) = Pn + l(z)+d n P n (z) _ {22Q) 

z 

For the recurrence coefficients in this case we have |27] 
I , b n +i — d n 



J n — "n 



b n — d n -i ' 



n = 1,2, ... 



d = d -oi, = d n _i ^ ra+1 ^" , n = l,2, (2.21) 

On — «n-l 

Note that there is some "irregularity" in the expression for d n in (|2.2ip for n = and n = 1, 2, . . . . 
This irregularity can be avoided if one formally puts 6o = 0. However in our case bo ^ 0. Hence 
we will indeed have such an irregularity in the analytic dependence of the coefficients d n in n. 
Namely, we have the explicit expressions 

■- - ( "+ 1/2)2 n=1.2 



(n+ l)(n + 2) ! 



) -"J ■ ■ ■ 3 



do = 1/8, d n = — , n = l,2,.... (2.22) 

n + 2 
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The corresponding moments are transformed simply as a "shift" 

cf ) = ^±i. (2.23) 

Comparing (|2.20p with (|2.16p we see that in our case the biorthogonal partners P n (z) coincide 
(CO 

with the CT LBP P„ for fj, = 0. Taking into account that in our case c± = — 1 we obtain that 
the "negative" moments are expressed as 

c_ n = -c n+ i, n = l,2,.... (2.24) 

The LBP are connected with the two-point Pade approximation problem [7]. Given the 
moments c n , n = 0, ±1, ±2, . . . , consider two formal power series 

oo oo 

F + (z) = ^c k z- k , F^{z) = Y J ^ k z k . (2.25) 

k=l k=0 

Then we have [27] 

^iM-^fl + OU— 1 ) ^M--EM + o{^) (2 26) 

P»(*) " c, +U[Z h P n (z) ~ Cl +U[Z j - (2 ' 2bj 

It is convenient to introduce the formal Laurent series 

= F + (z)-F_M (2 ^ 27) 

Cl 



Then we have 



" i^r " I o(zn, z - o, (2 - 28) 

i.e. LBP P n (z),P^ 1 (z) solve the problem of two-point Pade approximation (at z = 0, oo). Note 
that the Hermite formula (jl.lOp describes "one half" of this two-point Pade approximation 
problem. 

So it is reasonable to refer to the polynomials P n {z) as the Hermite elliptic Laurent biorthog- 
onal polynomials. 

3 The weight function and biorthogonality 

Consider the functions 

_ 2J(k) _ 2J(l/fc) 

F ' [z) -~WK{ky F+{z) -~lcoM> (3 - 1} 

where z = k 2 . It is assumed that the function F_(z) is defined near z = while the function 
-F+ (z) is defined near z = oo. Note that formally 

F + (,)-*ML. 

z 

From the considerations of the previous section it is easily verified that formula (|2.28p holds 

for the LBP defined by the recurrence coefficients (|2.3p . It is seen also that if one introduces 
the function (cf. [9]) 

_ F + (,)-F_(.-) (3 2) 

v ; 2mz v ' 
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then 

-k 



z~ k w{z)dz = c fc , fc = 0,±l,±2,..., (3.3) 



c 



where the integration contour C is the unit circle. Thus the biorthogonality property (|2.6p can 
be presented in the form 



P n (z)P m (l/ z)w(z)dz = h nm 5 nm . (3.4) 

c 

Now we calculate the weight function w(z) in a more explicit form. We have 

If J{k) J(l/k) \ 
W[Z) ~^{ ~k^K{k) + lWJk))' { } 

Using the relation J(k) = K(k) — E(k) and the formulas for the complete elliptic integrals with 
inverse modulus (as usual K'(k) = K(k'), E'(k) = E(k'), k' 2 = 1 - k 2 ): 

K(l/k) = k(K(k) + «.(*)), E Wk} = m-^-^KW + ^K'jk) 

k 

we arrive at the formula 

-K{k)K'(k) + E(k)K'{k) + K(k)E'(k) 
W[Z) ~ k*K(k)K(l/k) ' ( j 

The latter expression can be further simplified using the Legendre relation [23] 

-K(k)K'(k) + E(k)K'(k) + K{k)E\k) = tt/2 
to give 

W{Z) = *h K(k)K(l/kY (3 ' 7) 

Now we introduce variable 9 on the unit circle such that k = z 1 ! 2 = e 1 ®/ 2 . The biorthogonality 
relation can then be written as 

f27T 

P n (e ie )P m (e- ld )p(8)d8 = h n 5 nm , (3.8) 



where the weight function is 
i 

2e ie / 2 \K{e id l 2 )\ 

4 Generating function and explicit expression 



From (jl.lip we obtain the generating function for the corresponding LBP 

*(*>*) = m ^fe= m =^^1 (4.1) 
x 2 )(l- zx 2 ) ^ 



where ^ n is given by (12. 1]) and 

dt 



F(x;z) 



/0 y/(l-t 2 )(l-zt 2 ) 

is the standard (incomplete) elliptic integral of the first kind. 
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In order to find the explicit expression for the polynomials P n {z) from (|4.ip . we note that if 

^(i- X 2 ){1 - zx r ) n=Q 



X>n^ 2n (4.2) 



then, obviously, 

oo 

F(x;z) = Y,PnX 2n+1 /(2n + l) 



n=0 



and hence $>(x,z) = Yl A n (z)x 2n+1 where 

n=0 

Mz) = t^- (4-3) 

s=0 

Formula (14. 3D gives an explicit expression for the polynomials A n (z) and hence for the LBP 
P n (z) if the coefficients j3 n are known. But it is easy to verify (using e.g. the binomial theorem) 
that 

0n = i ^r L 2F 1 (-n, 1/2; 1/2 -n;z). (4.4) 
We thus have 



A (l/2) w (l/2) m (-n 
n{ ' ^(2m + l)n\m\(l/2- 



n (2m + l)n!m!(l/2-n) m 

m=0 

x 2F1 (-m, 1/2; 1/2 - m; z) 2 i ? i(m - n, 1/2; 1/2 + m - n; z). (4.5) 
Another explicit expression is obtained if one notices that 

1 f>V"y n (^-V (4.6) 



V(i-* 2 )U-^ 2 ) ^0 

where Y n (z) are the ordinary Legendre polynomials |15j : 

^n(t) = 2^1 f-n,n + l;l;^ 
We thus have the rather simple expression 

A n (z) = k"± Ys{ flf q) , (4.7) 

s=0 

where q = (k + k~ l )/2 (recall that k 2 = z). 

For z = 1 we have (3 n = \ for all ra in (|4.2|) . Hence we have from f|4.3|) 

A n (l) = G n , (4.8) 

where we denote 
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- the finite sum of inverse odd numbers. G n can be obviously expressed in terms of the Euler 
"harmonic numbers" defined as 



n 



H n = \^ l/k. 

k=l 

We have clearly 

G n = H 2n+1 - H n /2. (4.9) 

Consider the recurrence relation of the type (|1.5D 

(2n + l)ij) n - 2n(l + z)^n-i + (2n - l)^ n _ 2 = 0. (4.10) 

For z = 1 we see that ip n = G n is a solution of this equation. The second independent solution 
for z = 1 is trivial - it is a constant: ip n = const. This means that the general solution of the 
equation (|4,10p for z = 1 can be presented in the form 

^ n = a + PG n (4.11) 

with arbitrary constants a, (3. In particular, for B n (l) we can write 

B n (l) =G n -l. (4.12) 

5 Polynomials orthogonal on the unit circle 

Consider the Christoffel transform (|2.17p of our polynomials with fj, = 1. We have 

_ P w+1 (l) _ £ n A n+1 (l) _ (n + 3/2)G n+1 
n P n (l) Zn+iMV (n + 1)G„ " 1 ' ' 

For the corresponding transformed LBP, we have the expression 

PM = ^W-^-M, (5 .2) 

The moments are calculated by f|2. : 

Cra+l C n C n C n _|_i 

Cn " ci-1 " 2^ ' 



(5.3) 



Using property (|2.24p we see that c_ n = c n , i.e. the moments c n are symmetric with respect to 
reflection. In turn, this is equivalent to the statement that the corresponding polynomials P n {z) 
are the Szego polynomials which are orthogonal on the unit circle |21l [5]: 

2W P n (e ie )P m (e- ie )p(9)de = h n S nm , (5.4) 

o 

where 



P(0) = m r(e^-1). (5.5) 
ci — 1 



In our case we have explicitly (see (|3.9 
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It is well known that polynomials orthogonal on the unit circle are defined by the recurrence 
relation [S] 

P n+l {z) = zP n {z)-a n P*{z), 71 = 0,1,..., (5.7) 

where P^{z) = z n P n (l/z) (bar means complex conjugate). In our case all moments c n are real, 
hence -P^(z) = z n P n (l/z). The parameters a n are called the reflection parameters. They play 
a crucial role in the theory of Szego polynomials on the unit circle. We have 

a n = -P„+i(0). (5.8) 

It is well known [5] that if the reflection parameters are real and satisfy the condition \a n \ < 1 
for all n = 0, 1, . . . then the positive weight function p{9) > always exists. Moreover in this 
case the weight function is symmetric on the unit circle: p(2-7r — 9) = p{9). 
In our case we have 

a n = -P n+1 (0) = -U n+1 P n+1 (0) + P n+2 (0) = 1 - U n+1 , (5.9) 

where we used the property P n (0) = 1. 
It is easily seen that 

n = 0, 1,.... (5.10) 



2(n + 2) 2(ra + 2)G(n + 1) ' 

In (|5.10|) both terms are negative and less then 1/2 in absolute value. Hence — 1 < a n < for 
all n = 0, 1, . . . and the function p{9) is positive as is seen from (|5.6p . 

We thus have a (presumably) new example of the Szego polynomials orthogonal on the unit 
circle for which both weight function and recurrence coefficients are known explicitly. 

Following [4] and [26], to any polynomials P n {z) orthogonal on the unit circle with the 
property —1 < a n < 1 one can associate symmetric monic polynomials S n (x) = x n + 0(x n ~ 1 ) 
orthogonal on an interval of the real axis. Explicitly 

1 - a n -i 

where x = z 1 / 2 + z" 1 / 2 (it is assumed that one chooses one branch of the function z 1 / 2 such 
that for z = re %e we have z 1//2 = r 1 / 2 e lS / 2 > < 9 < ir). The polynomials S n (x) satisfy the 
three-term recurrence relation 

S n+ i + u n S n -i(x) = xS n (x), (5.12) 

where the recurrence coefficients are 

U n = (1 + On_i)(l - On_2), 71=1,2,.... (5.13) 

In (|5.13p it is assumed that a_i = — 1 (this is a standard convention in the theory of polynomials 
orthogonal on the unit circle), so U\ = 2(1 + ao). 

If the polynomials P n (z) are orthogonal on the unit circle 

i-2-k 

/ P n {e i6 )P^{e- td )p{9)d9 = 0, m + n (5.14) 
J o 

with the weight function p(Q) then polynomials S n (x) are orthogonal on the symmetric interval 
[-2,2] 

J S n (x)S m (x)w(x)dx = h n 5 nm (5.15) 
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with the weight function [H [26] 

where x = 2cos(#/2). 

In our case it is elementary verified that 

(n+1/2) 2 

u n = K 7 '{ , n = l,2,.... 5.17 
n{n + 1) 

These recurrence coefficients correspond to well known associated Legendre polynomials studied, 
e.g. in pp. Recall that generic associated Legendre polynomials are symmetric OP satisfying the 
recurrence relation (|5.12j) with the recurrence coefficients 

- = r IVm (5 ' 18) 

(n + vy — 1/4 

with arbitrary nonnegative parameter v. The ordinary Legendre polynomials correspond to 
v = 0. In our case we have v = 1/2. 

Consider the weight function for these polynomials. From (|5.16p and (|5.6p we derive 

W(X) = 2\K(e-^r (5 - 19) 

We can simplify this formula if one exploits the relations for K(z) where \z\ = 1. Namely, one 
has Q7] 

K{e ±i4> ) = ^e^ /2 {K{cos(<j)/2)) ± iK(sm((f>/2))), -vr < < vr, 
whence we have 

2 

W(X) = ^ 2 (cos(0/4)) + i^ 2 (sin(e/4))- (5 ' 20) 

Taking into account that x = 2 cos (9/2) we finally arrive at the formula 

2 

Ttf(ac) = ; ; -2 < X < 2. (5.21) 

iTVl/2 + x/4) + K*(y/l/2 - x/4) " " 

Note that the function is even w(—x) = w(x) as should be for symmetric polynomials. It 
has the only maximum at x = 0: w(0) = l//C 2 (y / l/2) = . Near the endpoints of the 



r 4 (i/4) ' 

interval [—2,2] the weight function w(x) tends to zero rapidly. 

The weight function for generic associated Legendre polynomials (with arbitrary v) was 
found in pp. Our formula (|5.2ip can be obtained from the results of pp by putting u = 1/2. 
Note, that Pollaczek studied [H] more general orthogonal polynomials containing 4 parameters. 
The associated Legendre polynomials (as well as the associated ultraspherical polynomials) are 
contained in the Pollaczek polynomial family as a special case. 

Similar polynomials were studied also in [19] where the author in fact rediscovered the Her- 
mite approach (as well Hermite's continued fraction (|1.9p ) to elliptic integrals. He consid- 
ered integrals J n (k) as moments for some "elliptic" orthogonal polynomials having w(x) = 
l/y (1 — x 2 )(l — k 2 x 2 ) as an orthogonality weight on the interval [—1,1]. As a by-product 
the author of [19] introduced polynomials which are similar to the Hermite polynomials A n (z), 
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B n {z). He then related them with the associated Legendre polynomials in a similar way. It is in- 
teresting to note that Hermite himself already introduced such polynomials in 111]. Hermite 
also established differential properties of these polynomials anticipating results of Rees |19j . 

Consider relations between LBP of special type and orthogonal polynomials on an interval 
in details. 

We have P n (z) = (P n+ \(z) — U n P n {z)) / '{z — 1). Substituting this formula to (I5.11|) and using 
the inversion property z n P n (l/z) = P n (z) of our LBP P n (z) we immediately obtain a very 
simple relation 

S n (x) =z~ n ' 2 P n {z). (5.22) 

We can see this also using recurrence relation for the LBP 

P n+1 (z) + d n P n (z) = z(P n (z) + b n P n ^{z)). (5.23) 

Assume (as in our case) that d n = —1, n = 0, 1, . . . . Then P n (0) = 1 for all n = 0, 1, 2, . . . and 
from (12. llh we obtain that the reciprocal polynomials coincide with the initial ones: 

z n P n (l/z)=P n (z). (5.24) 

Conversely, assume that some LBP are reciprocal (15.24j) . Then from (12. lip it follows that either 
d n = 1, n = 0, 1, . . . or d n = —1, n = 0, 1, . . . . But we have obviously P n (0) = 1 which leads to 
the only possibility d n = — 1, n = 0,1,2, ... . Thus condition d n = — 1 for all n is necessary and 
sufficient for LBP to be reciprocal invariant (|5.24p . In this case polynomials 

S n (x) = z~ n ' 2 P n {z) (5.25) 

are obviously monic polynomials in x = z 1 ! 2 + z" 1 / 2 . From recurrence relation (|5.23p with 
d„ = -1 we obtain recurrence relation for polynomials S n (x) 

S n+ i(x) + u n S n -i(x) = xS n (x), (5.26) 

where u n = —b n . Thus we arrived at the same symmetric polynomials on the interval as in 
the case of polynomials orthogonal on the unit circle. Note that for the LBP P n {z) with the 
property (|5.24p we have c\ = do = — 1 and c* = c n thus from (|2.12p we obtain 

Cl-n = "Cn (5.27) 

for all n = 0, ±1, ±2, .... Perform now the Christoffel transform with [i = 1 

P n (z) = P ^)-_UnPn(z) (5>28) 

From (|2.19p and (|5.27p we obtain that transformed moments are symmetric c_ n = c n . This 
means that polynomials P n (z) will satisfy the recurrence relation (|5.7p with reflection param- 
eters a n given by a n = 1 — U n+ \ = 1 — P n+ 2{1) / P n +\{1) ■ Formula (|5.1ip now is equivalent to 
formula (|5.25p . Thus starting from arbitrary LBP with the property d n = —1 we can arrive at 
the same symmetric OP on the interval S n (x) as for the case of polynomials orthogonal on the 
unit circle. From another point of view these relations are discussed also in [2]. 
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6 Geronimus transform. Laurent biorthogonal polynomials 
with a concentrated mass added to the measure 

In the previous section we showed that the Christoffel transformation of Hermite's elliptic LBP 
gives polynomials orthogonal on the unit circle with explicit reflection coefficients (|5.10j) and 
the weight function given by (|5.6p . In this section we consider another spectral transformation 
of Hermite's elliptic LBP which is called the Geronimus transform (GT). This transform was 
introduced in [27] and is similar to well known Geronimus transform for the ordinary orthogonal 
polynomials HUES]- Recall thatGT for LBP is defined as [27] 

P n (z) = V n P n {z) + z(l - V n )P n -i(z), (6.1) 

where 

V = l, V n = /- , n=l,2,.... (6.2) 

fi - (pn/Pn-l 



In (|6.2p /j, is an arbitrary parameter and <\> n is an arbitrary solution of the recurrence relation 

4>n+i + d n (j) n = n{(j> n + b n <t> n -i). (6.3) 

Note that (|6.3|) is the same recurrence relation that (|2.2p for LBP. Hence its the general solution 
(up to a common factor) can be presented in the form 

cj)n = Pn^)+XP { nlM- (6-4) 

It is easy to verify that the polynomials P n (z) are again LBP satisfying recurrence relation 

P n+1 (z) + d n P n {z) = z(P n (z) + ~b n P n ^(z)), (6.5) 

where the recurrence coefficients are 

h=x{Vi-l), b n = b n -i 1 — - , d n = d n ^p-. (6.6) 

1 — Vn-l Vn 



It can be shown that CT and GT are reciprocal to one another [27J. This observation allows 
one to obtain explicit biorthogonality relation for polynomials P n (z) starting from that for 
polynomials P n (z). Namely the pair of the Stieltjes functions F + (z), F—(z) is transformed 
as [27] 

~ uF + {z)+^ + u vF_{z) - jj-v 

F + (z) = — , F_(z) = — , (6.7) 

where 

MX MX 



ci-X d o-X 

Assume that < 1. Choose the contour C as the unit circle (if = 1 we choose C as the 
unit circle with a small deformation near z = [i in order to include the point z = \i inside the 
contour C). Then the weight function for new polynomials P n {z) is defined as in (|3.2[) . i.e. 

_ PM -P-W = «W + _%til_. (6. 8 ) 
Ziriz z — fj, 2itiz{z — fj,) 

The second term in (|6.8p will give a concentrated mass at the point z = jJL added to a "regular" 
part presented by the first term in (|6.8p . 
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Assume that [i = 1, as in our case. Then we have biorthogonality relation 

r-27T 



f p(e)Pn(e l °)P n (e- M )d9 + MP n (l)P n (l)=0, n + m, (6.9) 
J o 



where 



is the "regular" of the weight function on the unit circle and the last term ()6.10j) describes the 
concentrated mass 

M = ^±l (6.11) 

inserted at the point z = 1. 

In our case we have d n = —1, n = 0, 1, . . . hence v = — x/ (x + 1) an d 

M 



2jr(l + x) 



For the recurrence coefficients we have explicit formulas (16.60 where we need first to calculate 
the coefficients V n . We see that 

^n = J Pn(l) + X^ 1 _ ) l(l)- 

But the values P n (l) and P^\(l) were already calculated (see (14.80 . (14.120 ). We thus have 
(3/2)„ 



■(G n + 2x(G„-l)). (6.12) 



Now all coefficients V n , b n , d n are calculated explicitly. 

We thus constructed a nontrivial example of the Laurent biorthogonal polynomials with 
explicit both recurrence coefficients and the measure. The weight function for these polynomials 
has a concentrated mass on the unit circle. Note that in contrast to the Christoffel transformed 
Hermite's polynomials the polynomials P n {z) constructed in this section are not polynomials of 

the Szego type. This means, in particular, that the biorthogonal partners P n {z) do not coincide 
with P n (z). 



7 Associated families of the Laurent biorthogonal polynomials 

Return to the sequence J n (x,k) of incomplete elliptic integrals defined by (|1.2p . They satisfy 
three-term recurrence relation (jl.3p . Repeating previous considerations we can express J n (x] /c) 
for any n = 1,2, ... in terms of Jj(x;k) and Jj+i(x; k) for some fixed nonnegative integer j 
(in (jl.4p the case j = is chosen): 

J n+1 (x;k) = Q(x; n,j)y/(l - .x 2 )(l - Px 2 ) + A^Jj+^x; k) - B^J^k) (7.1) 

with some coefficients An \ Bn^ ■ 

We first note that there is an explicit expression of complete elliptic integrals J n (k) in terms 
of the Gauss hypergeometric function 

Jn(k) = 2^1(1/2, 1/2 + n; 1 + n; k 2 ) (7.2) 
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(relation (|7.2p can be easily verified by direct integration). 

Repeating similar considerations that were already exploited in the two first sections we can 
show that A$(z) and are determined by the recurrence relations 

iC?) iv„ i i\aU) „r„ i „■ i i/<r*/iO') 



(n + 3/2 + = (n + j + 1)(* + l)4f } - «(n + j + 1/2)4-1, 



(n + 3/2 + ijB^ = (n + j + l)(z + 1)B® - z(n + j + 1/2)^ (7.3) 



with the same initial conditions as (jl.6p . Hence, A^(z) polynomials of n-th degree in z = k 2 
and B#\z) is a polynomial of degree n — 1 multiplied by z. 
Introduce also the function 

where as usual we put fc 2 = z. 
Then we will have the property 

F(2;j) "SS =o(z " +i) - p - 5) 

n the same 



In the same way it is possible to derive the generating functions for polynomials An\z), 



k^J^x-k ) = y A U) {z)x 2n+l 

Using (|4.6p we arrive at the explicit representation for polynomials An (z) in terms of the 
Legendre polynomials Y n (q): 



Y s {q)Y n ^{q) 

s=0 

where z = k 2 , q = (k + l/fe)/2 



WW = (2j + y: vxrii » (7 - 7) 



Now we introduce the monic LBP P^z) = ^(z)/^ where £ n = (j + l)„/(j + 3/2) n . They 



satisfy the recurrence relation of type (|2.2p with 

<. = -!, t „ (" + W (7.8) 

(n + j)(n + j + 1) 

Thus we have j-associated polynomials with respect to the Hermite elliptic LBP, i.e. we should 
replace n — > n + j in formulas for recurrence coefficients. 

Due to the property d n = — 1 the associated polynomials P„ again possess the invariance 
property z n Pi j \l/z) = Pi j \z). We thus can construct polynomials orthogonal on the unit 
circle using formula (|5.28p . In order to get explicit expression for reflection parameters we 
need the value Pn\l). But this can easily be obtained from f|7.7j) : 

A^(l) = G n (j) (7.9) 
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where we introduced the function 

n 

(Up to a common factor G n (j) is a sum of n + 1 succeeding inverse odd numbers starting from 
l/(2j + 1); for j = it coincides with G n .) Thus 



P$i(l) n + 3/2+i G n (j) 



P#\l) n + j + 1 G n+1 (j) 



and 



£#> = 1 - 27, 



n+l- 



Corresponding symmetric polynomial Sn\x) on the interval satisfy recurrence relation (15, 26ft 
with 

(n + j + 1/2) 2 

n " (n + i)(n + ; + l) (7 ' iUj 

so they coincide with the associated Legendre polynomials considered in pQ. Indeed, from 

formula (|5.18p we see that the "shift" parameter v = j + 1/2 where j = 0, 1, 2, . . . . 

(i) 

It is interesting to note that recurrence relations for polynomials Pn can be presented in 
such a form that the all coefficients are linear in n - see e.g. ()7.3p . In [6] we considered a family 
of such LBP connected with so-called generalized eigenvalue problem on su(l, 1) Lie algebra. In 
our case, however, corresponding representations of su(l, 1) will not be unitary, in contrast to [6]. 
This leads to an interesting open problem how to describe associated classical LBP in terms of 
non-unitary representations of su(l, 1) algebra. Note that generic orthogonal polynomials with 
linear recurrence coefficients in n were studied in details by Pollaczek [18] who derived explicit 
expression for them and found the weight function as well. 

Note finally that LBP with recurrence coefficients (|7.8)1 belong to a family of so-called "asso- 
ciated Jacobi Laurent polynomials" introduced and studied by Hendriksen 0[9]. Nevertheless, 
in our case the parameters of the associated LBP belong to the exceptional class which was 
not considered in [51 [9] . This means that in some formulas in [8j [9] the bottom parameter in 
the Gauss hypergeometric function 2F1 (z) takes negative integer values. In this case formulas 
obtained by Hendriksen should be rederived in a different form. In particular, a simple explicit 
expression of the power coefficients (in terms of the hypergeometric function 4^3(1)) for as- 
sociated LBP obtained in [5J seems not to be valid in our case. Instead, we obtained explicit 
expressions like (|7.7p . 

8 Laurent biorthogonal polynomials connected 
with the Stieltjes— Carlitz elliptic polynomials 

Return to formula (jl.ip . If one denotes 

/•CO 

r„(p)= / k 2n sn 2n {t)e- pt dt (8.1) 
Jo 

then we obtain the recurrence relation 

p 2 r n = 2n(2n + l)r n+1 - 4(1 + k 2 )n 2 r n + 2n(2n - l)r n _i. (8.2) 
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Again it is seen that for every n > one can present 

r n+ i = A n n - B n r , (8.3) 

where obviously ro = e~ pt dt = p^ 1 and A n , B n satisfy the same recurrence relations as r n+ \ 
e.g. 

p 2 A n = 2(n + l)(2n + l)A: 2 A„_i - 4(1 + k 2 )(n + 1) 2 A„ + 2(n + l)(2n + 3)A n+1 (8.4) 

with initial conditions Aq = 1, A_i = 0, Bq = 0, -B_i = —1. Now it is seen that A n are 
polynomials of degree n in both variables p 2 and k 2 . If p = then ^4 n become polynomials in k 2 
introduced by Hermite. 

Coefficients A n considered as polynomials in p 2 become orthogonal polynomials because they 
satisfy three-term recurrence relation typical for orthogonal polynomials. Orthogonal polynomi- 
als of such type (and several related ones) where introduced and studied by Carlitz. He exploited 
some explicit continued fractions found by Stieltjes. These continued fractions are connected 
with elliptic functions (for details see, e.g. [13]). Today these orthogonal polynomials are known 
by Stieltjes-Carlitz elliptic polynomials |16^ I3"l [T3]. Note that the Stieltjes continued fraction is 
obtained from (|8.2|) by the same way as Hermite obtained his continued fraction (jl.9p for a ratio 
of two elliptic integrals. Consider now polynomials A n as LBP with respect to the argument 
z = k 2 . Passing from A n to monic polynomials P n (z) (by the same way as for the case p = 0) 
we arrive at the recurrence relation (12.21) with 



A 1 P 2 h in + l/2) 2 

dn -~ l ~W+W' K ~~ n(«+l) ' ( } 

We see that the recurrence coefficient b n is the same as for the Hermite LBP, but the coefficient d n 
now depends on n. This means that polynomials P n {z) do not possess symmetric property 
like (jZISj) . 

In contrast to the case p = the polynomials P n (z) have more complicated properties. For 
example, the reciprocal polynomials Pn(z) defined by (|2.10p do not belong to the same class, 
their recurrence coefficients appear to be 

- 1 (n+l)(n + 2) 2 (n + l/2) 2 

Cn ~ l + i ^' n " ((n + 2)2 +p74)((„ + l) 2 +P 2 /4)- 

The biorthogonal partners P n (z) defined by (12. 15ft have the recurrence coefficients 
n(n + l)(p 2 (n + 1) - n - 2) 



d n 

(n 4- 1 I^Vtr, 4- 1 \ 2 (r, 2 (n 4- 1 ^ - n - 0\ 

(8.7) 



{n{p 2 - 1) - l)((n + l) 2 + p 2 /4)((n + 2) 2 + p 2 /A) ' 

Hence biorthogonal partners P n (z) also do not belong to the same class that initial LBP P n {z). 

For the ordinary Stieltjes-Carlitz elliptic orthogonal polynomials the orthogonality measure 
can be found explicitly: it is purely discrete one (see [3l Q3]). The problem of finding the 
orthogonality measure for the corresponding LBP seems to be much more complicated. 
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